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Abstract 

New unstable temperature gradient driven modes in an inhomogeneous 
plasma are identified. These modes represent transient lo ~ f^W^th^^ sound 
oscillations in magnetized plasma that are kinetically destabilized via Landau 
interactions. Electron and ion sound branches are unstable for large values of 
the Larmor radius parameter k±pe,i ^ 1, respectively. The instability occurs 
due to a specific plasma response that significantly deviates from Boltzmann 
distribution in the region k±pi^e S> 1 • Pacs: 52.35 Kt, 52.35 Qz 
Small scale instabilities driven by temperature gradients are believed to be responsible for 
particle and energy transport in a tokamak. The ion temperature gradient driven turbulence 



1 



has been identified as a primary cause of ion tliermal transport. Tlie electron temperature 
gradient driven modes are believed to be largely responsible for the electron transport 0. 
Both types of modes, the ion temperature gradient (ITG) and electron temperature gradient 
(ETG), have been extensively studied over last decades In basic slab geometry these 
modes are essentially sound waves destabilized by coupling to pressure fluctuations 0-0- 
In toroidal geometry the modes are driven by the unfavourable magnetic curvature rather 
than by acoustic oscillations, though the mode may remain essentially slab-like in the region 
of the weak/negative shear [T^j2^. In this work we report on a new mechanism of the 



destabilization of the temperature gradient driven modes. The mechanism is essentially 
related to the wave-particle Landau interaction. However, the most important element is 
a newly found regime of plasma response for large values of the Larmor radius parameter, 
k±Pa ^ 1, « = In this regime plasma response in essential way deviates from usually 

assumed Boltzmann distribution for fc^p^ ^ 1. As is shown below modified plasma response 
occurs for transient (acoustic type) modes whose eigen- frequency u ~ k\\Vtha does not grow 
with increase of k^p"^, but rather remains constant. As a result two new branches of the 
temperature gradient driven modes occur for k\pj ^ 1 and k\pl ^ 1. We call them ion 
and electron kinetic temperature gradient (KTG) modes, respectively. 

We consider shearless dispersion equation in a local approximation. Within the local 
theory the parallel velocity and density perturbations for each species are given by standard 
expressions (see e.g. |[12| , p!7| ) 

SaDa, (1) 



Vtha T, 



a 



Here 

+ (1 + s^Ziso,)) (ro(6«) - ri(6,)) (3) 

UJ 



and 

L = i-(i-—) ro(&.) - ^ (ro(6.) - Ub^)) b^. (4) 

Various plasma parameters are defined as follows: = —kycTa/ CaB^Ln^ = 

-kycTa/eaBoLTa,L~^ = -riQ^duo/dx, Lxa = -T~^dTa/dx, s„ = uj/kuVtha, K = klpl/2, 
vf^^ = 2Ta/m„, Pa = Vtharriac/ (eaBo); To,i{b) = /o,iexp(-6), and Z{s) is the standard 
plasma dispersion function. We have introduced auxiliary potential = — uj/{k\\c)A, 
where is the electrostatic potential and A is the magnetic vector potential. 
The vector potential A can be found from the Ampere's law 

J|| = en{viii - ^e\\) = -^VlA. (5) 

By using (0) and (||) one obtains 

A- H{Der + Dd fc||C 

kl5^ - 2sl {DeT + A) a; ^' ^ ^ 



and 

2 _ . 



(7) 



for the parallel potential 0. Here r = T^/Ti, and 5^ = c^/Up^. 
By using the Poisson equation, 

- Vi0 = 47re(ni-ne), (8) 

and expressions (|^) and (|^ we obtain general local dispersion relation 

klS^ ikr + k + DiT + D,) - 2sl {D,T + D,) {kr + Q 
= -klXl {kl6' - 2sl {D,r + dS) , (9) 

where A|, is the Debye lenghth, A|, = Tg/ (47moe^). This general dispersion equation describes 
both ITG and ETG modes as well as new kinetic temperature gradient driven modes. 

To illustrate existence of a new unstable mode we plot a solution of the dispersion equa- 
tion @ as a function of the normalized ion Larmor radius parameter pj/L„. The following 



parameters are fixed kyPi = -\/2 x 0.3, k^pi = \f2 x 0.1, /? = 2 x 10~^, PilL^i = x 0.1, 
Pi/Lxe = X 0.01, fc||Pj = X 0.002, r = 1, We have chosen the above plasma param- 
eters to be exactly the same as in Ref. 17. The plasma pressure parameter j3 is defined as 
(3 = SirrioTi/BQ. [Note that our definition of fj^^, and (3 differ differ from those in Ref. 17 
by a factor of 2.] Figures la and lb show the real Ur and imaginary Ui parts of the eigen 
frequency, respectively. It is already evident from Figs la and lb that there are exist two 
distinct eigenmodes: in Fig. lb the left curve is the unstable ITG branch and the right 
curve is the unstable drift mode branch fl^ (in terminology of Ref. 17). Our results in Figs 
la and lb are in complete agreement with Ref. 17. 

To clarify the nature of the unstable branch we solve the dispersion equation as a function 
of the kyPi parameter for fixed (3 = 2 x 10'^, Pi/Ln = \pl x 10^^, Pi/LTi = \/2 x 10"^, 
Pi/Lj'f. = \pl X 10"^, k^pi = \pl X 10^^, r = 1. We also take A^i = for this case. The 
mode growth rate and mode frequency normalized to fcyfthe are shown in Figs. 2a and 2b, 
respectively. Two new unstable branches exist in the regions kyPi > 1 and kyPe > 1. The 
standard ITG mode has a peak growth rate around kyPi ~ 1, and the standard ETG mode 
has a growth rate peaked at kyPi ~ 40 (corresponding to kyPe ~ 1). A new ion mode (ion 
KTG) has a maximum growth rate around kyPi ~ 5, and a new electron mode (electron 
KTG) saturates to a constant growth rate quickly for large kyPi > 200 ( kyp^ > 5). It is 
essential to note that new modes exist in the regions of large Larmor radius parameters 
(respectively, kyPi for the ion KTG, and kyp^ for the electron ETG). As we show below, 
kinetic temperature gradient modes occur due to a peculiar behaviour of the plasma response 
in the region of large kyPa- Let us consider this response in more details. 

A standard notion is that for large values of the Larmor radius parameter, kyp^ > 1, 
the density response of the respective plasma component is Boltzmann due to decaying 
asymptotics of ro,i(&a) ~ l/y/K. for large ba, so that we have from (|]-|) Da <^ 1,Zq ~ 1. 
This, in fact, implicitly assumes that the ratio of ujna/^^ is finite for large kyPa- In turn, 
this requires that the mode eigen-frequency increases with kyPa (linearly or faster). The 
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latter is definitely true for drift wave type modes, where uj ~ uj^a- However, temperature 
gradient driven modes, are basically sound waves whose frequency is of the order of the 
transient (or sound) frequency, k\^Vtha- If the parameter k\^Vtha remains constant (so oj is 
approximately constant), the ratio oJna/^ will be increasing with kyPa that may compensate 
for the decaying factor from ro,i(6a). Then, for large kyPa ^ 1 we may write in the 

leading order 

-Toipo) = 7r^ — —r, (10) 



UJ 2^7? e LuLr, 



and 



^ (ro(&„) - TM) 6„ = -V-^- (11) 

By using these asymptotics for the electron KTG mode in the region kyPe ^ 1 and kyPi ^ 1 
we obtain the following response functions 

1 Vte 1 Vte /.„s 
le = l + (13) 



D, = LJ^I + ^J^I, (14) 

A^/tx ujLn sj 4v/vr ujLti sj ' 



+ J^-7^(l + ^e^(^e)). (15) 

The electron KTG mode is a destabilized electron sound wave with uj ~ k\\Vthe so we retain 
a full plasma dispersion function in De to account for the electron Landau interactions. In 
this region uj/k\\Vthi ^ 1 so we can simplify Z{si) functions. Solution of dispersion equation 
(y) with simplified response functions (P^P^) and A/? = is shown in Fig. 2 by squares. We 
note here, that in fact plasma response is not given exactly by expressions (111-0) but contains 



additional terms of the order of Pi/Ln |T^. These additional terms do not affect asymptotic 



expressions (|T^-|T^ as long as Pi/L^ < 1. 

A similar expansion can be made in the region kyPi ^ 1. Then the ion response for the 
ion KTG branch becomes 

k = l-^^ + ^^^ (16) 

A = irj=^ (1 + SiZ{s,i)) - ^ Si (lz{si) - Si - s^iZ{si)) 

^ [l + s,Zis,)). (17) 



2^/7^ uLTi 

The electron response functions can be simplified in this region by explicitly using k±pe < 
1 and uj/k\\Vte < 1- Then one obtains 

= 1 - (i - ^:^) , (18) 

D,= {l-^){l-2sl + ^s,V^) 1 + ^s, (-2.. + , (19) 

h = l-(l-^){l- K) - "^K. (20) 



D,= [l-^j[l-2si + ISeV^) (1 - K) + -^Se I -2Se + 1 " ^e) 



f^(l-2s2 + ZSeV^) fee. (21) 



Mode growth rate and real frequency obtained with (|l8]-0) and (^ with Ad = are 
shown in Fig. 2 by circles. 

Recently, it has been emphasised that the ETG modes can be significantly modified by 
a finite value of the Debye length parameter [0. Electron KTG mode is strongly stabilized 
in a high temperature plasma by the finite Debye length as shown in Fig. 3. Note that effect 
of the Debye length screening may be different for the negative shear plasma . 



Though the modified ion response for large ion Larmor radius is criticial for the ion KTG 
mode, the mode itself is destabilized by the electron temperature gradient as demonstrated 



in Fig 4. Ion KTG mode exists for large gradients of the electron temperarure (dashed 
line, r/e = 5). The mode growth rate is significantly reduced for smaller values of rj^ (solid 
line, rje = 0.5). Thus we may conclude that ion KTG mode is in fact a hybrid mode, in 
which both ion and electron response are essential. Electromagnetic effects are expected to 
be important for this mode existing in the intermediate region between standard ITG and 
ETG modes. Though for larger rje we observe a standard finite f3 suppression (curve labeled 
by diamonds in Fig. 5), for low rje = 0.5 , finite plasma pressure effectively leads to the 
extension of the ITG instability into the shorter wavelengths (curve labeled by squares in 
Fig. 5). It is interesting to note that for j3 = 0.1 taken for this case, the skin depth size 
c/upe ~ 0.1 Pi, which makes kyc/upe < 1 in the region of the unstable ion KTG. 

In summary, we have identified two new branches of the temperature gradient driven 
modes. These modes are closely related to acoustic type modes that are destabilized by 
the temperature gradient effects. A modified plasma response in the region of kyp^ ^ 1 is 
essential for new instability. There is certain similarity of the KTG modes to the universal 
instability of drift waves that is driven by density gradient and wave-particle interaction 
| |21| |. The KTG modes considered in this work will be affected by a finite value of the 
magnetic shear but noting a special role of the temperature gradients and the modified 
plasma responce for kyPa ^ 1 it remains unclear whether these modes can be stabilized by 



the shear similarly to the universal instability [p2| , p3| . One can expect that KTG modes can 
persist in the region of weak and/or negative magnetic shear where the mode structure is 
essentially that of a shearless slab [^. Toroidal effects will further modify the KTG modes 



but are not expected to cause their stabilization [24]. Generally, investigation of a general 



case requires nonlocal integral analysis [|T3|,|T5|,|T9[ that is left for future work. 
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FIGURES 

FIG.l. Normalized (to the ion cyclotron frequency) frequency (a) and growth rate (b) of 
the ITG and unstable drift wave modes versus the normalized ion Larmor radius pi/ L^. The 
solid hne corresponds to the ITG branch, while the dashed hne is that of the destabihzed drift 
wave branch. Parameters are the same as those in Ref. /3 = 2 x 10-^ pi/Lri = ^2 x 10-\ 
pi/LTe = V2 X 10-2, k^pi = x/2 X lO'S kyPi = 3\/2 x 10'^ k\\Pi = 2\p2 x lO'^. 

FIG. 2. Normalized wave frequency (a) and growth rate (b) for the ITG , ETG modes and 
kinetic modes. Plasma parameters: /3 = 2 x 10~^, PijLn — -\/2 x 10"^, pi/Lxi = x 10~^, 
Pi/Lxe = \/2 X 10"^ kxPi = \/2 X 10"^ k\\pi = 2\/2 x 10"^, r = 1. Sohd line - standard 
ITG and ion KTG; dashed hne - standard ETG and electron KTG. 

FIG. 3. Stabilization of the Kinetic Electron Temperature Gradient mode by the Debay 
screening effect: solid line - Xo/Pe = 0, dotted-dashed line - Xo/Pe = 0.70, dotted line - 
Xd/Pc = 2.21. Plasma parameters: /3 = 2 x 10"^, Pi/Ln = "\/2 x 10"^, Pi/Ln ^ V2 x 10^^, 
Pi/LTe = a/2 X 10-\ k^pi = a/2 X 10-1, A;||p, = 2a/2 x 10"^, r = 1. 

FIG. 4. Effect of the electron temperature gradient on KTG modes. Solid line - pi/ Ltc = 
\pl X 10"^ (77e = 0.5), dashed line Pi/Lxe — x 10"^ {rie — 5). Plasma parameters: 
/3 = 2 X 10-^ a/L„ = V2 X 10-2, p,/LTi = v^x 10-^ k^p^ = v^x 10-^ k\\pi = 2v^x 10-^ 
T = 1. 

FIG. 5. Effect of a finite plasma pressure on the KTG modes. Triangles - Pi/L^e — 
\/2 X 10-\ /3 = 2 X 10"^; diamonds squares - Pi/Lre = "\/2 x 10""^ (3 = 0.1; circles 

Pi/Lxe = V2 X 10-2, (3 = 2x 10-^ squares - pi/Lre = x 10-^, p = 0.1. Other plasma 
parameters: - pi/L„ = y/O.xlO''^, pi/Lxi = "\/2xl0-\ k^pi = "\/2xl0-\ k\\pi = 2"\/2xl0-^, 
T = 1. 
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Fig.2b A. I. Smolyakov 10/23/01 
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Fig.4 A. I. Smolyakov 10/02/01 
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